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Collisional  Relaxation  of  Transient  Vibrational  Energy  Distributions 
in  a Thermal  Unimolecular  System.  The  Variable  Encounter  Method. 

D.  F.  Kelley,  L.  Zalotaib  and  B.  S.  Rabinovitch 

Department  of  Chemistry  BG-10,  University  of  Washington 
Seattle,  WA  98195 

Abstract 

A method  has  been  developed,  called  the  Variable  Encounter  Method, 
for  the  study  of  the  relaxation  of  an  initial  vibrationally  cold  ensemble 
of  molecules  into  a vibrationally  hot  distribution  by  a known  and  variable 
number  of  successive  collisions  with  a hot  wall.  The  theory  of  the  experi- 
ment is  presented.  The  system  studied  was  the  isomerization  of  l,l,cyclo- 
propane-d^  with  a fused  quartz  wall  temperature  of  800  K to  1175K,  and 
average  number  of  collisions  from  2.3  to  22.3.  Various  modified  gaussian 
and  exponential  models  of  energy  transfer  were  found  to  give  agreement  with 
the  data.  The  average  down-step  size  was  found  to  decline  from  $ 3500  cm’^ 
at  the  lowest  temperature  to  % 2500  cm”^  at  the  highest  on  the  basis  of  a 
gaussian  model.  A mathematical  analysis  of  the  relation  between  mean  first 
passage  times  and  incubation  times  is  given.  Incubation  times  increase 
from  ^ 7 to  ^ 12  collisions  with  increasing  temperature.  Transient  population 
distributions  and  the  sequential  reaction  probabilities  as  a function  of 
collision  number  are  calculated. 
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I.  INTRODUCTION 

There  have  been  no  previous  measurements  of  the  dynamics  of  the  collisional 
relaxation  of  vibrationally  excited  polyatomic  molecules.  If  an  ensemble  of 
reactant  molecules  embedded  in  a hot  bath  gas  starts  out  in  a "cold"  distribu- 
tion, whether  Boltzmann  or  not,  it  will  relax  by  collisions  with  the  bath  gas 
molecules  to  a steady-state  "hot"  distribution.  The  study  of  reaction  proba- 
bilities associated  with  the  intermediate  distributions  can  provide  insights 
into  the  form  of  the  energy  transfer  distribution  function  over  the  entire 
range  of  energies  from  0 to  Eq,  the  reaction  threshold.  One  may  thus  determine 
the  entire  matrix  of  collision  energy  transfer  probabilities  P . The  proba- 
bility matrix  also  appears  in  expressions  for  the  mean  first  passage  time, 

tfp,  of  molecules  with  respect  to  an  absorbing  barrier,  as  given  by  Kim1  and 
2 

Widom.  In  the  trivial  case  of  a constant  temperature,  steady-state,  low 
pressure,  unimolecular  system,  t^p  = l/kQ,  where  kQ  is  the  low  pressure  rate 
constant. 

Formal  theoretical  treatments  of  the  transient  in  vibrational  relaxation 

3 4 

have  been  given  for  special  assumed  forms  of  P.  ’ Usual  treatments  have  con- 
sidered  the  relaxation  to  be  a continuous  rather  than  discrete  process.  This 
is  valid  for  the  case  of  a fixed  time  interval  between  collisions  only  if  the 
relaxation  time  is  long  compared  to  the  time  between  collisions.  This  condition 
is  not  met  for  a reasonably  strong  collider.  A general  treatment  then  requires 
that  the  relaxation  process  be  formulated  in  terms  of  difference  rather  than 
differential  equations.  The  resultant  difference  equations,  which  are  not  tract- 
able by  any  analytical  technique,  must  be  solved  numerically. 

A very  simple  method  of  study  of  the  collision-by-collision  relaxation  of 

5 6 

vibrational  energy,  i.e.,  of  the  relaxation  transients,  has  been  developed. 

The  collisions  are  with  a wall  which,  in  analogy  with  experimental  results  on 
homologous  series  of  bath  molecules  of  increasing  size,^  and  on  the  basis  of  a 
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quasi -stati stical  model,  is  anticipated  to  be  at  least  as  efficient  as  any  gas- 
eous bath  molecule.  In  this  technique  a molecule  experiences  a large  number  of 
collisions  with  a cold  wall  of  known,  variable  temperature  and  is  brought  into 
equilibrium  with  it  before  having  a series  of  collisions  (an  encounter)  with 
the  hot  walls  of  a reactor.  By  varying  the  geometry  of  the  reactor,  one  can 
change  the  average  number  (and  number  distribution)  of  successive  hot  collisions, 
hence  the  name.  Variable  Encounter  Method  (VEM).  Escape  from  the  reactor,  usually 
without  reaction,  leads  to  re-equilibration  on  the  cold  wall  before  re-entry 
into  the  hot  reactor.  The  amount  of  reaction  that  occurs  after  a given  number 
of  encounters  may  be  measured. 

The  reaction  chosen  for  this  study  is  the  ring  opening  isomerization 

of  1 ,l-cyclopropane-d2  to  propene.  This  process  has  a critical  threshold  of 

15  -1 

about  64  kcal  and  a pre-exponential  factor  of  ^ 1.3x10  sec  . It  has  been 
studied  in  the  homogeneous  steady  state  case,10  and  a preliminary  account  of  VEM 
work  has  been  reported.®  Light  cyclopropane  has  also  been  used  in  VEM  work  to  be 

5 

reported.  This  system  was  chosen  not  only  because  homogeneous  studies  of 
energy  transfer  at  comparable  temperatures  with  inert  bath  gases  already 
exist,  10  but  also  because  this  reaction  is  well  known  to  be  notably  free 
from  complicating  surface  and  wall  seasoning  effects. 
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II.  EXPERIMENTAL 


The  apparatus  consisted  of  a 5-liter  fused  quartz  spherical  flask  which 
had  several  quartz  cylindrical  finger  reactors  of  various  dimensions  blown  onto 
the  surface.  Each  finger,  which  could  be  heated  separately,  consisted  of  a 
cylindrical  region  of  length  £ , terminated  by  a hemisphere  of  radius  r 
which  varied  from  2.1  to  2.5  cm  for  different  reactors.  The  ratios  of 
£/r  for  the  reactors  used  are  given  in  Table  I.  The  entrance  area  to  a 
reactor  was  •v  1%  of  the  total  area  of  the  flask.  The  reactor  could  be  enclosed 
in  a stainless  steel  heater  block  which  was  heated  by  clamshell  Kanthal  ceramic 
heaters.  Fine  silica  sand  was  used  as  a thermal  conduction  medium  in  the  approxi- 
mately 2-3  mm  gap  between  the  quartz  reactor  and  the  steel  block.  Temperature 
measurements  were  made  with  several  chromel-alumel  thermocouples  glued  onto  the 
outside  of  the  reactor  finger.  Typical  reactor  temperatures  were  750K  to  1150K, 
while  flask  temperatures  were  350K  to  400K. 

Some  temperature  gradients  over  the  reactor  surface,  particularly  near 
the  attachment  to  the  cold  flask,  were  unavoidable.  This  problem  was  sub- 
stantially reduced  by  having  additional  heaters  embedded  in  the  top  of  steel 
block.  In  this  way,  the  reactor  temperature  could  be  maintained  constant  to 
± 5°  within  a distance  of  'v  0.7  cm  from  the  top  of  the  reactor  where  the 
temperature  conmenced  a more  abrupt  decline.  The  temperature  profile  over 
the  longer  reactors  (where  radiation  losses  were  not  so  great)  could  be  main- 
tained constant  within  2°-3°.  Thus,  an  uncertain  area  of  reduced  temperature 
amounted  to  n*  10-15%  for  the  smallest  reactor  and  was  negligible  for  the  longer 
reactors.  Calculations  indicate  that  temperature  gradients  through  the  quartz 
wall  2 mm  thick)  were  completely  negligible. 


A pre-exposure  of  the  oxidized  quartz  surface  to  the  substrate,  at  a pressure 
-4 

of  about  5x10  torr  deposited  a thin  film  of  "polycyclopropane"  on  the  surface 
and  provided  reproducible  results  for  the  measured  rates  of  isomerization.  The 
seasoning  time  required  to  achieve  reproducible  results  was  somewhat  longer  than 
in  conventional  static  systems  (i.e.,  hours  and  days,  rather  than  minutes  or 

hours).  Prior  to  a run,  the  system  was  pumped  down  to  < 10”6  torr,  the 

-4 

reactant  was  introduced  at  a pressure  usually  equal  to  1-2x10  torr,  although 

-4 

experiments  were  made  over  the  range  0. 5-9.0  x 10  torr.  Run  times  varied 
from  ^ 10  minutes  to  several  hours  with  the  amount  of  reaction  being  1-50%. 

Analysis  was  by  gas  chromatography  with  a squalane  SCOT  column  and  a 
flame  ionization  detector. 
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111 . CALCULATIONAL  TECHNIQUE 

1 . Calculation  of  “number  of  collisions"  distribution  per  encounter 
There  is  a distribution  n of  the  number  of  collisions  made  by  a 
molecule  with  the  wall  after  entering  the  reactor  (an  encounter).  The  distri- 
bution n was  changed  by  changing  the  fc/r  ra-io;  in  this  work,  £/r  varied  from 
0.17  to  10.  This  distribution  was  obtained  by  a stochastic  calculation  of  the 
acutal  trajectories  of  single  molecules.  In  the  calculation  of  n , it  was 
assumed  that  the  surface  is  rough  on  the  molecular  level,  so  that  a cosine  law 
distribution  of  angular  reflection  is  valid  after  each  collision  regardless  of 
incident  trajectory.  (If  this  assumption  were  replaced  by  one  of  specular 
reflection,  it  would  result  in  only  a modest  change  in  n . The  details  of  this 
calculation  are  now  given. 

Molecules  which  cross  the  boundary  plane  separating  the  reactor  finger  from 
the  bulb  do  so  with  an  angular  distribution  function 

P (0 ) = 2 sine  cose  , (1 ) 

where  P(6)  is  defined  as  the  normalized  probability  of  a molecule  trajectory 
making  an  angle  6 with  respect  to  the  normal  to  the  plane.  This  cosine  law 
distribution  is  the  general  result  for  effusive  molecules.  Actually,  at  typical 
operating  pressure,  the  mean  free  path  is  comparable  to  the  dimensions  of  the 
spherical  flask.  However,  the  same  result  is  obtained  by  integrating  over  the 
surface  of  the  flask  (see  App.  A). 

The  z axis  is  chosen  along  the  axis  of  the  cylinder;  z = 0 defines  the 
plane  separating  the  hemispherical  and  cylindrical  parts  of  the  reactor  finger, 
with  z positive  in  the  hemisphere.  The  xQ  and  yQ  coordinates  on  the  zQ  = -i.fr 
plane  (separating  the  cylinder  from  the  large  flask)  are  chosen  randomly  for 
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any  one  molecule,  such  that  xQ  + yQ  5 1 , in  units  of  r.  The  initial  trajec- 
tory of  the  molecule  is  defined  by  xQ  and  yQ,  by  6 , chosen  randomly  to  conform 
to  eq.  1,  and  by  an  azimuthal  angle  <p  , chosen  randomly  between  0 and  2it.  This 
trajectory  may  be  written  as 


(x(t) ,y(t) ,z(t) ) = (x0.y0.z0)  + 


where  t is  a (continuous)  parameter,  and  x‘ = sine  cos<j> , y'  = sin6  sin$, 
z'  = cose.  This  trajectory  may  intersect  the  reactor  in  either  the  cylinder 
or  the  hemisphere;  the  two  cases  are  handled  separately.  In  either  case,  one 
calculates  the  first  impact  coordinates  (x^,y^,z^).  To  obtain  the  new  set  of 
coordinates  (x",y",z"),  which  define  the  next  trajectory  analogous  to  eq.  (2), 
one  first  defines  a transformation  such  that. 


0 « T y 


the  new  coordinates  are  then  given  by. 


1-COS0  \ 
sine  simp  J , 
sine  cos <p/ 


where  6 and  $ are  chosen  in  the  same  manner  as  for  the  Initial  trajectory. 

For  the  case  of  the  impact  on  the  hemisphere,  the  matrix  representation  of 

T is  given  by  a product  of  rotations, 

(cosip  0 sinA/cosa  sina  o\ 

0 1 0 Jl-sina  cosa  0 1 ; 

-simp  0 cosip/\  0 0 1 / 

ip  = sin  1 (z^;  a = tan  ^(y^/x^.  For  impact  on  the  cylinder,  T is  a product 
of  a rotation  about  the  z axis  and  a translation.  The  next  impact  point  is 
found,  with  consideration,  also,  of  the  possibility  of  escape  from  the  finger. 
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This  process  is  repeated  until  the  molecule  escapes.  In  this  way,  the  distri- 
bution n is  determined  by  calculating  many  trajectories  (1000-5000,  practise). 
An  average  inter-collision  free  flight  distance  is  also  obtained. 

From  n , it  is  useful  to  define  a correlation  table  for  each  reactor. 


i=j+l 


C.  is  defined  as  the  probability  that  a molecule  will  have  at  least  j+1  colli- 
J 

sions  given  that  it  has  j,  i.e.,  the  probability  that  a remaining  molecule  will 
not  escape  after  the  jth  collision  (Table  I).  The  average  number  of  collisions  per 
encounter,  m,  is  recovered  from  the  correlation  table  by  the  expression: 


m = l 
j = l 


After  a few  collisions,  the  molecule  "forgets"  the  initial  conditions  and  the 
correlation  table  asymptotically  approaches  a constant.  The  number  of  colli- 
sions required  to  reach  this  limit,  within  statistical  fluctuations  introduced 
by  the  stochastic  nature  of  the  calculation,  varies  from  0,  for  the  shortest 
reactor,  to  about  15,  for  the  longest. 


2.  Simulation  of  the  "collision"  process.  Models  for  energy  transfer 
A normalized  Boltzmann  population  vector  N , characteristic  of  the 
flask  temperature,  is  constructed  to  represent  the  molecules  that  enter  the 
reactor.  To  simulate  a collision,  this  vector  is  operated  on  by  the  probability 
matrix.  After  each  wall  collision,  the  population  of  each  molecular  energy  level 
changes  until  an  equilibrium  distribution  (or  steady  state  if  reaction  is 
considered)  is  realized.  The  population  of  the  ith  level  after  collision, 

N?  , is  given  by 


Ni 


l p..N. 
• 111 
J=0  J J 
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where  N.  is  the  population  of  the  jth  level  before  the  collision,  and  p..  is 
J * j 

the  probability  per  collision  of  a transition  from  the  jth  level  to  the  ith 
level.  This  equation  may  be  written  in  matrix  form  N'  = P N . Since  no 
change  takes  place  after  an  equilibrium  distribution  of  molecules  is  established, 

j I (P  “ 1)  Neq  = 0 (3) 

The  following  constraints  thus  apply  to  the  P matrix, 

Pi j/Pji  = ^9i/gj )exp[-(Ei~Ej )/RT] 

00 

1 p. . = 1 

. it 

where  gi  is  the  degeneracy  of  the  ith  level. 

Since  there  are,  of  course,  far  too  many  levels  to  permit  their  enumera- 
tion, the  levels  are  grouped  together,  usually  in  "grains"  of  100  cm  . 

If  we  now  specify  the  form  (model)  of  the  down  transitions,  this  will  completely 
specify  the  matrix  P . Further  details  of  the  construction  of  this  matrix  are 

*w> 

given  in  ref.  (11).  In  the  calculations  two  different  models  for  the  proba- 
bility of  a given  down  transition,  AE  , were  used: 

exponential:  P^^  = A exp(-AE/<AE>)  ; 0 < AE  s 9000  cm-1  ; 

=0  ; AE  > 9000  cm'1  ; 

gaussian:  P^E  = A ' exp[- (AE-AEmp)2/2o2]  ; 0 5 AE  5 9000  cm'1  ; 

=0  ; AE  > 9000  ; 

where  the  limitation  AE  < 9000  cm’1  is  a practical  computational  feature  to  limit 
the  total  size  of  the  matrix  to  tractable  dimensions;  A and  A*  are  normalization 
constants  which  also  depend  on  the  up  transitions;  <AE>,  AEmp  (mp  signifies  most 
probable),  and  o are  parameters  of  the  model.  We  note  that  on  the  exponential 
model,  if  <AE>  « 9000  then  the  effective  size  of  the  average  down  transition  <AE> 

L k ...  ... — 
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is  simply  <AE>.  For  the  gaussian  model,  if  o « AEmp  then  this  effectively  becomes 
a stepl adder  model  and  <AE'>  = AEmp;  if  o l AEmp  then  there  is  a fairly  large 
probability  for  elastic  or  near-elastic  transitions,  and  this  form  begins  to  resemble 
a poisson  model. 

In  the  case  of  the  exponential  model,  <AE>  (or  AEmp  in  the  case  of  the  gaussian) 
was  also  allowed  to  vary  with  initial  energy  in  some  calculations.  This  dependence 
was  given  by  <AE>  = a + 8EX,  with  the  same  expression  for  AEmp,  where  a,  8,  and  x 
are  adjustable  parameters.  Most  thoroughly  studied  were  the  cases  where  8 = 0,  a 
adjustable  (a  "flat"  model),  and  a = 500,  8 adjustable,  x = 1 (a  "linear"  model). 

(a  is  not  chosen  to  be  zero  in  the  linear  model  to  avoid  a singularity  in  the  con- 
struction of  the  matrix).  In  all  cases,  down  transitions  larger  than  the  initial 
energy  were  treated  as  elastic.  The  probability  matrices  so  constructed  conform  to 
eqs.  3 and  4. 

In  the  limit  of  a strong  collider,  one  expects  a Boltzmann  distribution  to 
result  from  every  energy  level  of  the  initial  distribution, in  one  collision.  To 
transform  the  molecules  initially  at  high  energy  into  a Boltzmann  distribution 
requires  a larger  down  step  than  is  required  for  molecules  initially  at  lower 
energy.  Thus, linear  models  have  some  pretension  to  greater  physical  content 
than  do  flat  models. 


3.  Simulation  of  the  entire  encounter  process 

Between  successive  wall  collisions,  say  the  jth  and  j+lth,  the  population 
vector  must  be  attenuated  for  loss  due  to  reaction  and  escape.  The  ith  element 
in  the  vector  is  modified  by 

Ni  = cj e Ni 

where  is  the  jth  element  in  the  correlation  table;  ki  is  the  microscopic  RRKM 
rate  constant  at  energy,  Ei  ; and  t is  the  average  intercollision  time  obtained 
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from  the  average  intercollision  distance,  assuming  translational  velocity 
characteristic  of  surface  temperature.  The  vibrational  parameters  used  in 
the  RRKM  calculation  of  ki  as  a function  of  energy  are  given  in  App.  B. 

This  collision-attenuation  process  is  repeated  until  the  vector  N.  e 0 for 
all  i,  or  until  some  extrapolation  procedure  may  be  used  (see  App.  C).  The 
amount  of  reaction  after  each  collision  is  then 

-k.t  -k.x' 

P.1  * l N.CCjd-e  1 )+  (1  - CjHl-e  1 )] 


where  t‘  is  the  average  flight  time  before  a gas-gas  or  cold  wall  collision 
after  a molecule  leaves  the  reactor  finger,  both  treated  on  a fixed  gap  basis. 
In  this  system,  the  effective  pressure  is  sufficiently  low  so  as  to  make  the 
calculation  of  P'  insensitive  to  the  fixed  gap  approximation,  or  to  any 

J 

reasonable  error  in  the  calculation  of  the  k^.  The  average  probability  of 


Figure  1 presents  the  experimental  average  probability  of  reaction 
per  collision  Pc (m ) versus  temperature  for  the  five  reactors.  Least-squares 
quadratic  fits  are  shown;  the  corresponding  Arrhenius  plots  are  given  in  Fig  2. 
Values  of  the  observed  activation  energies  are  summarized  in  Table  II.  If  the 
wall  were  a strong  collider,  there  would  be  no  transient,  and  all  Pc (m ) curves 
would  coincide.  (Related  behavior  in  strong  collider  low-pressure  thermal  uni- 
molecular  theory  is  that  the  rate  of  reaction  of  substrate  embedded  in  a bath 
gas  at  temperature  T is  independent  of  the  initial  energy  distribution  of  sub- 
strate molecules  below  the  critical  threshold  EQ.)  Figure  3 shows  the  least 
squares  fit  curves  and  comparison  with  several  calculated  curves.  In  all  cases, 
<AE 1 > was  adjusted  to  fit  the  m = 2.3  curve,  at  temperatures  where  such  data 
was  available. 

Several  features  of  the  calculated  curves  reveal  themselves: 
a)  The  calculated  value  of  Pc(m)  is  a more  sensitive  function  of  <AE’>  for 
the  smaller  reactors.  This  is  the  reason  why  <AE ' > was  in  all  cases  adjusted 
to  fit  the  data  from  the  m = 2.3  reactor  where  available.  (It  should  be  kept 
in  mind,  however,  that  this  reactor  is  probably  the  one  most  susceptible  to 
experimental  error,  i.e.,  temperature  measurement  and  gradients  over  the  surface.) 
b)  For  small  m,  i.e.,  the  case  for  the  first  few  collisions  when  relatively 
more  of  the  molecules  are  at  low  energy,  the  exponential  model  is  more  efficient. 
This  is  due  to  its  relatively  greater  probability  of  very  large  up  jumps;  this, 
in  turn,  permits  a larger  contribution  to  reaction  after  only  a few  collisions 
than  do  the  gaussian  distributions  so  that  a relatively  smaller  value  of  <AE> 
is  required  in  order  to  fit  the  m = 2.3  data.  For  large  m,  by  contrast,  more 
molecules  have  moved  up  to  higher  energies  near  E (Figs.  4-6)  and  gaussian 
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models  give  larger  values  of  F (m)  than  do  exponential  models;  the  reason  for 
this  is  the  high  probability  of  near-elastic  collisions  associated  with  the 
exponential  model.  Similar  effects  may  also  be  seen  (Fig.  3)  upon  comparison 
of  the  two  gaussian  models  characterized  by  different  o,  where  the  one  with 
larger  a has  properties  more  closely  related  to  an  exponential  distribution. 

Linear  models  give  higher  values  of  Pc(m)  than 
do  the  corresponding  flat  models  for  large  m reactors.  The  reason  for  this  is 
that  the  linear  model  has  a lower  <AE*>  at  low  energies  (Figs.  4-6)  and  a smaller 
probability  of  effecting  transitions  from  these  to  energies  above  EQ  (so  that  a 
larger  value  of  <AE‘>  is  required,  to  fit  the  data  for 

the  m = 2.3  reactor).  For  large  m,  the  higher  <AE‘>  associated  with  the  linear 
model  at  higher  energies  enhances  the  reaction  probability  compared  to  the  flat 
model.  The  two  effects  mentioned  above  produce  the  result  that  the  linear 
exponential  and  flat  gaussian  (o  = 0.7  AEmp)  curves  are  virtually  coincidental. 

As  can  be  seen  from  Figure  3,  three  different  models  all  give  adequate 
(within  experimental  error)  fit  to  the  data.  The  best  fit  seems  to  be  obtained 
with  the  flat  gaussian  (o  = 0.7  AEmp)  model  and  the  linear  exponential  model. 

The  flat  gaussian  models  with  o = 0.35  AE^  give  poorer  fit  and  are  somewhat 
above  the  experimental  curves.  The  linear  gaussian  models  (not  shown)  lie 
slightly  above  the  corresponding  flat  models,  especially  at  higher  temperature, 
c)  Regardless  of  the  model,  two  findings  emerge:  i)  <AE‘>  drops  significantly 
with  increasing  temperature;  and  ii)  <AE’>  for  wall  collisions  is  greater 
than  <AE>  for  gas-gas  collisions  (see  Table  III).  This  drop  of  <AE>  with 
temperature  is  the  reason  why  the  Pc (m ) curves  separate  more  at  higher  tempera- 
tures. At  lower  temperatures, one  loses  resolution  between  the  models  as  they  all 
approach  closer  to  the  strong  collider  behavior. 

i A 
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2.  Sequential  reaction  probabilities  and  population  distributions 
The  form  of  Pc(m)  may  be  written  as 

Fc<m>  - s jl,  P<J)  ci  o-p<n) 

where  is  the  ith  entry  in  the  correlation  table;  thus  the  continued  product 
represents  the  fraction  of  reactant  molecules  left  after  j-1  collisions;  P(j) 
is  the  sequential  reaction  probability,  i.e.,  the  probability  of 
reaction  after  the  jth  wall  collision  made  by  a molecule.  Plots  of  P(j)  for 
various  models  are  given  in  Fig.  7.  The  smaller  m reactors  weight  the  small 
j values  of theP(j)  more  than  do  the  larger  m reactors.  Thus,  additional  expla- 
nations to  the  effects  discussed  in  the  previous  section  may  be  seen  in  terms 
of  these  plots.  The  reason  that  the  m = 2.3  reactor  is  much  more  sensitive  to 
changes  in  <AE‘>  than  the  larger  in  reactors  is  due  to  the  fact  that  P(j)  is 
much  more  sensitive  to  changes  in  <AE‘>  at  small  j than  at  large  j.  That  is, 
the  rate  of  relaxation  is  much  more  sensitive  to  changes  in  <AE>  than  is  the 
value  to  which  it  relaxes.  This  has  been  found  to  be  true  for  all  models. 

The  relative  and  absolute  population  distributions , which  formed  some 

basis  for  the  discussion  in  the  previous  section,  are  shown  in  Figs.  4-6. 

Examination  of  these  figures  shows  that  the  population  distribution  is  very 

significantly  altered  after  one  collision.  As  a consequence,  for  these  fairly 
strong  collisions,  the  reaction  rate  should  be  quasi -independent  of  input 

distribution  (flask  temperature).  This  was  found  to  be  the  case  experimentally 

when  the  flask  temperature  was  raised  by  100°  without  effect. 

It  is  also  clear  from  these  distributions  that  the  average  energy  approaches 

its  steady  state  value  rather  rapidly.  This  is  due,  in  part,  to  the  fact  that 

there  is  a very  large  average  upstep  at  low  energies.  The  large  upstep  is  a 

result  of  detailed  balance  at  low  energy  levels  where  the  densities  of  molecular 
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states  change  very  rapidly  with  energy.  A related  conclusion  is  that  the 

"temperature",  as  determined  by  the  bulk  of  the  population  distribution, 

relaxes  much  more  quickly  than  does  the  population  around  EQ,  i.e.,  than  does 

the  reaction  rate.  In  the  case  of  the  gaussian  model  (o  = 0.7  AEmp)  at  high 

temperature,  this  upstep  is  so  large  that  the  value  of  Emp  after  one  collision 

overshoots  the  steady  state  value  a little  (Fig.  4).  Despite  the  fact  that 

this  gaussian  model  works  well,  we  question  how  physically  reasonable  this 

result  can  be.  Indeed, as  o was  made  smaller,  other  non-physical  cal culational 

results  were  obtained  at  lower  energies  with  the  gaussian  models.  For 

o = 0.17  AEmp,  P(j)  ceased  to  be  a monotonic  function  of  j.  For  the  present 
system,  in  which  the  energy  range  of  interest  varies  from  zero  upward,  a simple 
step  ladder  model  is  quite  unsatisfactory.  This  is  in  contrast  to  chemical 
activation  systems,  where  only  energies  around,  or  above,  EQ  play  an  important 
role.  A variety  of  additional  population  distributions  for  all  models  discussed 
have  been  calculated  and  may  be  found  in  ref.  6b. 

The  relative  decline  in  the  Arrhenius  activation  energy  with  decrease  in 
m (Table  II)  illustrates  the  relative  effect  of  increasing  population  depletion, 
at  higher  energies,  with  rise  of  temperature  in  the  transient  regime. 

3.  Mean  first  passage  time  and  the  "incubation"  time 

.1  2 

The  first  mean  passage  time,  t^p,  as  developed  by  Kim  and  Widom 
gives  the  average  amount  of  time  required  for  a molecule  to  attain  an  energy 


equal  to  or  greater  than  Eq  as 


f°° 

t p(t)  dt  ; 
Jn 


y\or  nn  a.  I i x e<L 

where  p(t)  is  the^probabi 1 ity  that  a molecule  will  cross  Eq,  for  the  first 
time,  at  time  t.  t^p  is  a molecular  property  and  is  independent  of  m. 
Clearly,  tfp  is  a well-defined  quantity  at  all  pressures  but  is  most  simply 
related  to  the  actual  rate  of  reaction  in  the  low  pressure  limit  of  a uni- 
molecular  system.  If  the  system  is  at  low  pressure,  p(t)  is  the  probability 
of  reaction  at  time  t,  which  is  the  case  we  now  consider  and  to  which  our 
system  approximates. 
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Equation  5 may  be  written  as  a sum 

°°  j-1 

• l jP(j)  n (i-P(i)) 

Tp  j=l  i=o 

where  time  is  measured  in  units  of  collision^,  and  the  continued  product  is 
the  fraction  of  unreacted  molecules  after  j-1  collisions,  At  t s u.m(y 
We  define  an  "incubation"  time  in  terms  of  mean  first  passage  times, 

r tSS 

Tinc  ~ rfp  " rfp 

where  tmsfs  is  the  steady  state  mean  first  passage  time,  and  tfp  is  the  mean 
first  passage  time  for  the  case  of  transients.  It  is  shown  in  Appendix  0 that 

Tinc  = 1im  (j  ' P(i)/P)  + °(j+2p)  (6) 

j -*•  °°  i=o 

where  P is  P(°°),  and  j+  is  a characteristic  relaxation  time  for  the  system, 
and  P(j+)  > 0.6P.  In  most  cases,  the  second  term  of  eq.  6 is  small;  in 
the  present  case,  it  is  n.  10"^  - 10"^  of  the  first.  So  for  this 
case,  we  may  write 

j-1 

tI  = lim  (j  - l (P(i  )/P) ) (7) 

j -*■ 00  i=o 

Dove  and  Troe  ^ define  an  incubation  time  in  shock  systems  by  the 

equation 

N(t)/N(o)  = exp[-P(»)(t-T^c)]  (8) 


where  N(t)  is  the  concentration  of  reactant  at  time  t,  and  where  eq.  (8) 
is  to  be  evaluated  at  sufficiently  large  t so  that  all  transients  have 
disappeared.  It  is  shown  in  Appendix  D that  t^c  is  simply  the  special 
case  Tjnc  (eq.  7).  Calculated  values  for  xjnc  are  given  in  Table  IV. 
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Conclusions 

For  the  temperatures  studied  in  this  system,  relaxation  approaches  the 
steady  state  in  a small  number  of  collisons  'v  10-20  (Fig.  4)  with  correspond- 
ing incubation  times  of  7-13  collisions  (Table  IV).  This  accounts  for  the 
finding  of  steady  state  behavior  in  low  pressure  pyrolysis  studies1^  where 
m 5 80-100. 

Adequate  agreement  with  the  data  is  obtained  with  flat  exponential, 
linear  exponential,  and  flat  gaussian  (o  =0-7AEmp)  models.  The  best  agree- 
ment is  with  linear  exponential  and  flat  gaussian  models  whose  calculated 

P (m)  curves  are  essentially  indistinguishable.  The  flat  gaussian  (o  = 0.35AE  ) 

c mp 

linear  gaussian  and  stepladder  models  do  not  give  adequate  agreement. 

Wall  collisions  are  more  efficient  than  gas  phase  collisions;  a quartz 
surface  seasoned  with  a film  of  "polycyclopropane"  is  characterized  by  down 
step  sizes  of  ^ 3800  and  2950  cm’1  (flat  gaussian  model.  Table  III)  at  800K 
and  975K,  respectively.  These  compare  with  lesser  values  of  ^ 3300  and 
1900  cm  1 deduced  previously  on  a stepladder  model  for  neat  cyclopropane^ 
gas  phase  collisions  at  the  same  temperatures  (note  these  values  would  decline 
somewhat  on  a gaussian  model  basis). 

The  data  of  this  study  also  support  the  previous  finding1^  that,  in  the 
temperature  range  investigated,  energy  transfer  collisional  efficiency  declines 
at  higher  temperatures  in  polyatomic  molecule  systems. 
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Table  I.  Correlation  Tables,  C , for  Various  Reactors 

n 


l/r 

0.17 

0.46 

1.45 

4.0 

10.0 

n\n 

2.3 

2.9 

5.0 

10.5 

22.3 

1 

0.566 

0.659 

0.750 

0.755 

0.759 

2 ' 

.566 

.660 

.788 

.829 

.833 

3 

.661 

.800 

.882 

.889 

4 

.661 

.809 

.901 

.908 

5 

.815 

.915 

.922 

6 

.819 

.923 

.933 

7 

.822 

.928 

.942 

8 

.824 

.930 

.950 

9 

.824 

.931 

.958 

10 

.931 

.963 

11 

.968 

12 

.972 

13 

.975 

14 

.978 

15 

.979 

16 

.980 

17 

* 

1 
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Table  II.  Apparent  Arrhenius  Activation  Energies 

m 

■ ' — ■ i . ■ — . . . .1-  ..  . I". ' ■ 

Li  LI  M 10.5  22.3 

( kca 1 mole  ^ ) not  defined  35.5  42.5  48  50 


Table 

III. 

AE  Values  for 

Various 

Models 

(cm-1 ) 

T(K) 

Model 

8251 

900 

950 

1025 

1100 

1150 

flat  expon- 

<AE> 

VI 0000 

4300 

4000 

3000 

2550 

2300 

ential 

<AE’> 

^ 6500 

3035 

2940 

2685 

2280 

2115 

linear  expo- 
nential2 

<ae> 

^22000 

6900 

5890 

4765 

3870 

3420 

<ae'> 

^ 7800 

3500 

3400 

3160 

2900 

2720 

6 

% 1.0 

0.285 

0.24 

0.19 

0.15 

0.13 

flat  gaussian 
o=  0.35  AE 

AEmp  ! 

3500 

3200 

3100 

2930 

2700 

2530 

mp 

<AE'>  ! 

flat  gaussian 

AE 

3200 

2800 

2700 

2540 

2330 

2175 

o=  0.7  AE 

mp 

mp 

<AE ' > 

3500 

3100 

2990 

2810 

2580 

2410 

1 AE's  at  this  temperature  were  selected  to  fit  m = 10.5  curve  since  data  for 
m = 2.3  was  not  available.  Also,  for  exponential  calculations,  the  maximum 
AE  (matrix  truncation)  is  chosen  to  be  18,000  cm"^  and  the  grain  size  to  be 
200  cm”^;  otherwise  the  customary  9,000  cm"^  and  100  cm  ^ were  used. 

2 

AE  values  given  at  the  level  of  energy,  Eq  . 


18 


Table  IV.  Incubation  times,  t,  (number  of  collisions) 


*nc 

Model 

900K 

1150K 

flat  gaussian  (o=0.74E  ) 

mp' 

6.5 

10.7 

flat  exponential 

7.4 

12.3 

1 inear  exponential 

8.0 

13.3 

I 
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Appendix  A.  Derivation  of  Angular  Distribution  Function  for  Entry  in  the  Reactor 


Consider  a spherical  polar  coordinate  system  (r,0,$)  centered  in  the 
large  flask.  The  surface  separating  the  reactor  finger  from  the  flask  (which 
for  this  purpose  is  taken  to  be  a differentially  small  plane  of  area  A)  is  at 
the  point  0 = 0;  the  radius  is  normalized  to  unity.  We  examine  the  low  pressure 
case  where  gas-gas  collisions  may  be  ignored  and  derive  the  angular  distribution 
of  molecules  hitting  A assuming  a cosine  law  distribution  of  molecules  coming 
off  the  surface  of  the  flask. 

Now  consider  molecules  that  leave  the  sphere  from  the  angular  ring  0 to 
0 + d0  and  hit  A.  This  ring  has  area  2irsin0  d0.  The  angle  between  A and  the 
origin  of  the  coordinate  system  as  seen  from  any  point  on  the  ring  at  0 is 


given  by 


= (tt-0  )/2 


The  fraction  of  molecules  that  leave  this  ring  and  hit  A is 


df  = 2cost|i  dft/2ir 


where  dft  is  the  solid  angle  subtended  by  A.  (Note  that  if  dft  is  set  equal 

rV2 

to  the  "volume  element"  2irsini^  d^  then  df  = 1 so  eq.  A2  is  normalized.) 

0 2 

For  the  geometry  under  consideration,  dft  = Acost^/r  , where  r is  the 
distance  between  the  ring  and  A,  and  r = sin0/sin\j>  , so  that 

2 

A 2,  sin  \p 
df  = - cos  ip  ■ y 

* sin  0 

This  fraction,  multiplied  by  the  area  of  the  ring  and  divided  by  A,  gives 
the  fraction  of  molecules  incident  between  ip  and  i|/  + dip , 

P(*)d*  = I (2,  Sine  § d*) 


V s™  0 
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Using  eq.  A1 , this  simplified  to, 

P(t^)dijj  = 2 cosiJj  s i nv|>  dip 

which  is  the  cosine  distribution  with  the  appropriate  volume  element,  as 
desired.  Although  this  distribution  was  derived  assuming  A to  be  infinitesma! , 
it  holds  approximately  for  finite  A,  provided  A is  sufficiently  small  so  that 
the  element  of  spherical  surface  of  the  flask  may  be  approximated  as  a plane 
on  A. 
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Appendix  B.  Vibrational  Frequencies  for  RRKM  Calculations  (cm"1). 


1 ,l-Cyclopropane-d2  molecule: 

3100,  3080,  3040,  3020,  2330,  2210,  1480,  1350,  1180,  1130,  1110  (2) 
1070,  .1030,  1020,  980,  850,  810,  770,  620,  590. 
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Appendix  C.  Extrapolation  Method 


The  total  probability  of  reaction  per  encounter  may  be  written  as 

rj-l 


Penc  = mPc(m)  = I P(i)  I .n  ci(1-p(i)) 


j=l 


(Cl) 


1=0 


where  C.  is  the  ith  element  in  the  correlation  table  (C  =1); 
j-1  1 ° 

JI  C.(l-P(i))  is  the  fraction  of  molecules  that  have  at  least  j collisions. 
i=o  1 

P(i)  is  the  sequential  reaction  probability. 

P(j)  is  a monotonically  increasing  function  of  j that  asymptotically 
approaches  a limiting  value,  P(<»).  For  large  enough  j,  called  j*,  the 
approach  to  P(°°)  can  be  described  by  the  difference  equation. 


P ( k+ j * ) 


P(k+j*-l ) 


(C2) 


where  R .*  and  R .*  . 
J J ^ 


are  defined  by  the  following  equations 


P(j*-1)  = P(j*-2)  (1  + Rj*_1 ) , 


and,  similarly, 

P(j*)  = P(j*-l)(l  + Rj*) 


One  does  the  iterative  simulation  up  to  j*,  calculates  R.*_,  and  R.*  , and 

then  uses  the  extrapolated  values  of  P(j)  in  Cl.  It  is  found  that  eq.  C2 

workfquite  well  if  one  chooses  j*  such  that  R.*  5 0.01.  This  corresponds 

J 

in  the  present  study  to  a j*  of  about  15  to  30.  We  note  also  that  eq.  C2. 
implies  a simple  exponential  relaxation  to  P («>)  when  R^+j*  <<c  1 • 
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endix  D.  Incubation  Time  From  First  Mean  Passage  Time. 


The  incubation  time  is  defined  as  the  difference  of  mean  first  passage 
times  with  and  without  the  transient. 


t * t tSS 
Tinc  tp  xfp  * 


For  time  measured  in  units  of  collisions. 


l jP(j)  n (l-P(i))  - l jP(l-P)j  , 
j=l  i=o  j=l 


where  P is  P(°°). 

These  sums  may  be  performed  directly;  however,  they  converge  very  slowly 
so  that  this  becomes  computationally  undesirable.  The  following  technique 
may  be  used  to  approximate  the  sum  to  any  desired  level  of  accuracy.  We 
define  j‘  such  that  P(j')  has  virtually  converged  on  P and,  recognizing  that 
the  second  term  contains  a binomial  expansion  of  1/P  , we  rewrite  eq.  D1  as 

j ' j-1  oo  ••ij'-l 

•tinr  « l jP(J)  n (l-P(i))  + P l j(l-P)J'J  n (l-P(i))  - (l-P)/P.  (D2 ) 

1nc  j=l  i=o  j=j'+l  i =o 


We  note  that  t^  = ( 1 -P )/P  = 1/P  (where  t^  is  not  exactly  1/P  because  an 
integral  was  replaced  by  a discrete  sum).  Adding  and  subtracting  to  the  second 
term  of  eq.  D2, 


y J-1  j'-l 

Tinr  = l jp(j)  n (1-P(1))+Pn  (1-P(i))(l-P)-J  ( 

1nc  j=l  i=o  i=o 

L_ 


i))(l-P)'j’  (1-P)  l j(l-P)j_1 
j=l 


(1-P)/P  . 


j=l  J 

j-1  / 

(f-1 

) n (l-P(i ))- P 

n 

i=o  \ 

i=° 

j'-l 

n 

L 

i=o 

(l-P(D)  - (l-PJ/P 
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Of  course,  eq.  D3  must  be  exactly  equal  to  eq.  01  If  one  lets  j*  -*•  • 

(so  P(j')  = P exactly),  but  is  calculationally  tractable  because  the  sum 

converges  for  relatively  small  values  of  j';(in  the  first  term  of  eq.  D3, 

terms  in  the  sum  approach  zero,  and  in  the  second  term,  factors  in  the  product 

5 ft 

approach  unity,  as  j'  ■*■<*>•,  evaluation  of  eq.  D1  would  require  10  to  10  terms 
whereas  eq.  D3  requires  only  30  to  60).  Equation  D3  may  be  simplified  further 
by  expanding  in  a power  series  in  P (all  P(j)  are  considered  to  be  of  order  P). 
After  some  manipulation  with  use  of  the  binomial  expansion  of  (1-P)  ’J  and 
the  relation 

j‘-l  - j'-l  j'-l  i-1 

n d-P(i))  = 1-1  P(i ) + l P(i ) l P(i ' ) + •••  . 
i=o  i=o  i=o  i-o 


eq.  D3  becomes 


T . 

inc 


1 im 

. I 


j'-l 

j'  - l P(D/P 
i=o 


+ R + 0(P4) 


(D4) 


where 


j'  j’-1 

R=  lim  R..  = lim  £ j (P( j )-P)  + j ' ( j ' -1 )P/2  + (1-j ' ) l P(i) 
j ' ■*•*>  J j j=l  i=o 

j'-l  1-1 

+ l P(D  ,1  P(1  )/P  . 

1«o  i*o 

2 

Clearly,  R is  of  order  P.  If  P were  large,  such  that  the  0(P  ) term  should 
be  evaluated,  then  eq.  03  should  be  used.  If  P Is  sufficiently  small, 
eq.  04  becomes  simply 


Ttnc  ■ },'•  - f'  PtO/P) 

j -*®  i=o 


(05) 


-5  -8 

In  the  present  study,  P is  10  to  10  , so,  as  will  most  often  be  the  case, 


eq.  05  is  entirely  adequate. 
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2 

To  examine  the  intermediate  case,  we  ignore  all  terms  0(P  ) and  derive 

oo 

a more  explicit  estimate  of  the  magnitude  of  R.  We  write  R « l (R.-R  ); 

j=l  J J_l 

in  terms  of  this  sum  R is, 

J*  . j-2 

R.=  1 ini  l [j(P(j)-P(j-lJ  + P(j-l)  - P - (1  -P(j-1)/P)  l P(i )]  (06) 

j ' j=l  i=o 

R may  be  readily  calculated  for  the  important  special  case  of  exponential 

relaxation  of  the  transient,  i.e.,  P(j)  = P(l-e~^+).  If  the  relaxation  is 

not  too  fast,  then  we  may  replace  the  sums  in  eq.  D6  by  Integrals.  (Replacing 

the  sum  by  an  integral  in  eq.  05  in  this  case  gives  j+  = T'jnc-)  The  first  two 

00  ( J ' 

terms  of  eq.  06  cancel  and  we  are  left  with  R - P( j ' )/P-l ) P(j)djdj‘,  or 

R - - i j+  P.  In  general,  we  find  R is  0(j+  0),  where  j'  is  the  character- 
istic relaxation  time  of  the  transient.  The  exact  relationship  between 
and  the  j'  necessary  for  convergence  of  eq.  03  varies  with  the  nature  of  the 
relaxation  process. 

For  the  definition  of  given  by  eq.  8,  it  follows  that  (with  time 

measured  in  collisions,  see  (Fig.  8), 

Tinr  = lim  j’+p-ln(N(j,)/N(0) 
inc  j'--  P 

But 

j'-1 

- lnN( j )/N(0)  = l P(j)  ; 

j=o 

therefore, 

,DI  . ]1«  fj'  - Jf'  P(.i)/p|  ■ t'  , 
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Figure  Captions 

Fig.  1 Plots  of  experimental  values  of  Pc(m)  vs  T (K)  for  each  of  the 

five  reactors.  Solid  curves  are  quadratic  least  square  fits  to 
each  set  of  data.  The  m-number  is  attached  to  each  set. 

Fig.  2 Observed  rate  constants  calculated  from  least  square  fit  curves  of 

Fig.  1 plotted  y£  reciprocal  of  absolute  temperature.  The  m = 2.3 
curve  has  too  much  curvature  for  the  apparent  activation  energy  to 
be  well  defined. 

Fig.  3 Comparison  of  the  experimental  Pc(m)  vs^  T least  square  fit  curve 

(solid)  for  each  reactor  (Fig.  1)  with  curves  calculated  on  the  basis 
of  various  models:  — • — , flat  gaussian  (o  = 0.35  AEmp); 

, linear  exponential  and  flat  gaussian  (o  = 0.7  AEmp) 

, flat  exponential.  The  linear  exponential  and  flat  gaussian 

(a  = 0.7  AEmp)  curves  are  essentially  indistinguishable.  The  m = 10.5 
curves  are  not  shown  to  avoid  cluttering  of  the  graph;  however,  the 
relative  positions  of  these  curves  are  very  similar  to  those  for  the 
m = 22.3  case.  For  each  model,  the  AE  value  was  chosen  at  a particular 
temperature  so  as  to  fit  the  m = 2.3  curve;  the  values  are  given  in 
Table  III. 

Fig.  4 Histograms  of  the  calculated  sequential  reaction  probability  P(n)  ys^  n, 
the  number  of  consecutive  collisions.  Calculations  performed  at 
a)  900K  and  b)  1150K  with:  ...  , gaussian  (o  = 0.7  AEmp); 

, linear  exponential;  and  , flat  exponential  models, 

with  AE  values  indicated  in  Table  III. 
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Fig.  5 a)  Absolute  population  distributions  N(E)  vs  E at  900K,  portrayed 
after  various  numbers  of  collisions  (attached  to  each  curve)  and 
calculated  with  linear  exponential  or  (indistinguishable)  flat 
gaussian  (o  =0.7  AEmp)  models.  The  curves  on  the  right  represent 
magnification  of  ordinate  scale  by  a factor  of  100.  The  n = 17  dis- 
tribution closely  approximates  the  steady  state  distribution, 
b)  Absolute  population  distributions  N(E)  vs^  E at  1150K  calculated 
with  the  flat  gaussian  (o  = 0.7  AEmp)  model  and  portrayed  after  various 
numbers  of  collisions.  The  n = 21  distribution  curve  is  essentially 
steady  state.  See  caption  (a). 

Fig.  6 Absolute  population  distribution  at  1150K  calculated  with  the  linear 
exponential  model  after  various  numbers  of  collisions.  The  n = 25 
distribution  approximates  the  steady  state.  See  caption  Fig.  5a. 

Fig.  7 Relative  population  distributions  N(E)/N(E)gq  vs  E(cm~^)  calculated 

with  the  linear  exponential  model  at  a)  900K  and  b)  1150K,  portrayed 
after  various  numbers  of  collisions  (attached  to  each  curve).  In  both 
cases,  the  distribution  of  the  last  collision  shown  approximates  the 
steady  state  distribution. 


Fig.  8 Graphical  representation  of  eq.  8.  Solid  curve  is  asymptotic  to 
dashed  line  whose  slope  is  steady  state  reaction  rate  P . 


